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Abstract. We present a brief overview of a splinter session on determining the 
metallicity of low-mass dwarfs that was organized as part of the Cool Stars 16 confer- 
ence. We review contemporary spectroscopic and photometric techniques for estimat- 
ing metallicity in low-mass dwarfs and discuss the importance of measuring accurate 
metallicities for studies of Galactic and chemical evolution using subdwarfs, creating 
metallicity benchmarks for brown dwarfs, and searching for extrasolar planets that are 
orbiting around low-mass dwarfs. In addition, we present the current understanding of 
the effects of metallicity on stellar evolution and atmosphere models and discuss some 
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of the limitations that are important to consider when comparing theoretical models to 
data. 



1. Introduction 

Low-mass dwarfs are the most numerous stellar constituents of the Milky Way and 
have main sequence lifetimes that exceed the current age of the Universe (at least for 
those that are not brown dwarfs). They therefore form an important laboratory for 
probing the structure and evolution of the Milky Way's disks. Because of their ubiq- 
uity, cool dwarfs may represent the largest population of stars with orbiting planets, 
especially low-mass planets in their respective habitable zones, which are considerably 
closer for cool dwarf systems. In addition, the diminutive sizes of these stars makes the 
detection of transiting planets easier than for higher mass stars (for any given planetary 
radius). Previous results hav e demonstrated that pla nets are more likely to be found 



orbiting metal-rich stars (e.g. iFischer & Valentill2005l) . There were preli minary indica 



tions that the M dwarf s with known planets had sub-solar metallicities (IBonfils et al. 



|2b05; Bea n et all l2006h . in stark contrast to their high-mass counterparts. However, re 



cent results have shown that the M dwarfs with attending planets appear to be metal-rich 
(see Section 2; Johnso n & Apps 2009i). With low-rtiass dwarfs becoming imp ortant 



sites for planet hunting (e.g. MEarth; Irwin et al]|2009l : lEndl et aL ||2003; Johnson et all 



l2007b . the observational efficiency of these searches could be vastly increased with 
prior knowledge of stellar metallicity. 

Because the ages of low-mass dwarfs span the lifetime of the Milky Way, they 
can provide important insight into the history and evolution of the Galaxy. Recent 
improvements in kinematic modeling and magnetic activity analysis have provided en- 
hanced statistical age estimates for populations of low-mass dwarfs (West et al. 200i, 



20081) . Coupled with metallicity information, these ages can provide valuable insight 



into the chemical evolution history of the Milky Way disks. Without large samples of 
low-mass dwarfs, the utility of the statistically derived ages is limited. Fortunately, 
the advent of large surveys such as SDSS and 2MASS has pro duced photometric sam 



pies of low-mass dwarfs that number in the tens of millions (Bochanski et al. 2010) 
and spectroscopic samples that con tain more than 70,000 M dwarfs (West et al.,.2008l : 
Kruse et al.ll2010l : IWest et al.ll2010h and ahnost 500 L dwarfs (ISchmidt et al.ll2010bh . In 



addition, these large catalogs of low-mass dwarfs have identified significant samples 
of metal-poor subdwarfs. The detailed metallicities of these objects, coupled with their 
kinematic distributions, establish important constraints on the structure and composi- 
tion of the Milky Way halo. 

Historically, the metallicity of low-mass dwarfs has been an elusive fundamental 
property due to the complex atmospheres of M, L and T dwarfs that have restricted 
the accuracy of detailed model atmospheres. Over the past several years, new observa- 
tional techniques as well as independent theoretical advancements in atmospheric mod- 
els have produced results that appear to link the metalli city of low-mass d warfs to both 
their photometric and spectroscopic properties (e.^.. Bean et al. 2006; Bonfils et alj 



2005'; IWoolf & WallersteinI l2006l : iJohnson & AppsI 12009. : .Hauschildt & Baron. .201 



Rojas-Ayala et al. 2010). While these relations and the resulting metallicities provide 
fundamental measurements for stellar astrophysics, they also play a crucial role in stud- 
ies of Galactic evolution and the environments that host extrasolar planets. 
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Figure 1. The low-mass target stars of the California Planet Survey in the Mk 
\s V - K plane (filled circles), and the M dwarfs known to harbor one or more gas 
giant planets (five-point stars). The isometallicity contours for [Fe/H] = (solid) 
and [Fe/H] - +0.2 (dotted line) are ba sed on the broad-band photometric metalUcity 
calibration of iJohnson & AppsI (|2009|) . 



2. Calibrating M Dwarf Metallicity using Photometry 



Several previous studies have estimated M dwarf metallic ities using wide binary pairs 
that consist of both an M dwarf and a higher mass star (e.g., Bea n et al .12006': 'Bonfils et all 
2005 ; Woolf & Wallerste in 2006). Because binaries are assumed to be both coeval and 
have the same metallicity, the composition of the higher mass star (which can be accu- 
rately derived from comparison to theoretical models) can be applied to the companion 
M dwarf. Some of these studies have used opti cal and infrared s pectroscopy to tie 
spectroscopic features t o a metallicity scale (e.g.. iBean et al.ll2006l : IWoolf et al.ll2009l : 
Rojas-Ayala et alJbolOl) . Although the spectroscopic method has shown great promise 



in deriving M dwarf metallicities (see Section 3), spectroscopy is considerably more 
time consuming than photometry and may not be easy to obtain for large samples of 
starSj 

Bonfils et al. (iooi) used M dwarfs in wide binaries to derive a relation between 



the absolute A'-band magnitude and the V-K color (higher metallicity M dwarfs are 
slightly brighter at a given color). Given the large number of M dwarfs for which there 
exist photometric observations, this relation may prove exceedingl y usefu l. However, 
there were 2 problems with resulting analyses: 1) using the Bonfil s et al. 1 dloos) rela- 
tion, planet hosting M dwarfs appeared to be metal poor compared to their FGK star 
counterparts; and 2) and the relation yielded a mean metallicity of M dwarfs in the so- 
lar neighborhood that was almost 0.1 dex below the me a n [Fe/ HI of higher mass stars. 
These discrepancies were resolved by Ijohnson & Avvs ('2009'), wh o discovered a sys- 
tematic uncertainty in the photometry used by Bonfils et al.i (i2005l) . IJohnson & AppsI 
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Figure 2. A linear combination of the EWs of the Ca I and Na I features versus 
the H2O-K index for northern 8 pc M-dwarfs. The black dots represent M dwarfs 
with photometric metallicities and the yellow dots represent M dwarfs with only 
near-infrared (NIR) spectroscopic metallicities. The big black dots (with circles) 
represent the M dwarf planet hosts. Typical errors in EWs and H2O-K index are 
represented by the error bars. The dashed lines in the top panel are iso-metallicity 
contours for [Fe/H] values of -0.30, -0.05 and -h0.20, calculated from the NIR [Fe/H] 
calibration. The NIR [Fe/H] calibration allows to cover a larger sample of cooler and 
distant M dwarfs (yellow dots) 



(120091) used corrected photometry to r e-derive a relation between th e metallicity, Mk 
and V - K color of M dwarfs (see also lSchlaufman & Laughlinll2010l ). 

Figure[T]shows the low-mass target stars of the California Planet Survey in the Mk 
vsV-K plane (filled circles), and the M dwarfs known to harbor one or more gas giant 
planets (five-point stars). The isometallicity contours for [Fe/H] = (solid) and [Fe/H] 
= +0.2 (dotted line) are based on the broad-band photometric metallicity calibration 
of I Johnson & A pps (2009). The distribution of stars illustrates the tendency of planet- 
hosting M dwarfs to be metal-rich compared to stars in the Solar Neighborhood. The 
planet-metallicity relationship therefore holds for M dwarfs as well as Sun-like FGK 
stars. 



3. Calibrating M Dwarf Metallicity using Infrared Spectroscopy 



Most of the attempts to estimate the overall metal content of M dwa r fs have been per- 
forme d at visible wavelengths (e.g.. lGizislll997l : [Bonfils et alJl2005l : Ijohnson & AppsI 
2009h . Since M dwarfs are optically faint, this limited past analyses to early-type 



M dwarfs and few specific near by stars, which are brig ht and have accurate paral- 
laxes. To avoid this limitation, iRojas-Ayala et all (I2OIOI) developed a near-infrared 
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(NIR) [Fe/H] spectros copic calibration u sing strong absorption features in the ^T-band 
spectra of M dwarfs. Roias-Avala et al. (2010) adopted a similar approach to that of 



Bonfils et all (|2005) and Johnson & Apps (2009), assuming that binary systems share 



the same metallicity since both components formed from the same original molecular 
cloud. Seventeen FGK+M binary systems in the SPOCS catalog dValenti & Fischer] 
12005) were used as metallicity calibrators. The NIR [Fe/H] calibration uses the Equiv- 
alent Widths (EWs) of the Na I doublet and the Ca I triplet, and a water absorption 
index (H7O. ICovey et al.ll2010l) to differentiate between metal-rich and metal-poor M 
dwarfs (cr~0.15 dex). The results obtained with the NIR spectroscopic [Fe/H] cali- 
bration are in agreemen t with the results obtained with the photome tric calibration by 
iJohnson & Apps' ("2009"). The eight M dwarf planet hosts analyzed by Ro jas-Ayala et alj 
(12010) have metallicities higher than -0.05 dex, with the Jovian pl anets hosts being 
more metal-rich that their Neptune analogs. This corroborates the iJohnson & AppsI 
(l2009h conclusion that planets are found preferentially around metal-rich stars, like in 
their Sun-like counterparts. 

As a moderate resolution iiT-band spectrum can be efficiently obtained for most M 
dwarfs with current spectrographs (e.g. TripleSpec, FIRE), the NIR [Fe/H] calibration 
allows observations of cooler and distant M dwarfs (Figure |2]). Thus, this technique 
will enable the identification of likely planet hosts at lower masses than is possible with 
optical [Fe/H] techniques. However, the NIR [Fe/H] calibration is currently limited to 
M dwarf spectral types earlier than ~ M7 and [Fe/H] > -0.7 , due to the lack of FGK 
dwarf /late-type M dwarf wide binary systems with measured spectroscopic metallici- 
ties, and subdwarfs with A/AAk3000 NIR spectra, to be used as calibrators. 



4. Spectral Features of Low-Metallicity Brown Dwarfs 

Brown dwarfs are expected to have a similar metallicity distribution to the stellar com- 
ponents of our Galaxy, but reliably determining the chemical compositions of individ- 
ual brown dwarfs is a difficult task. Atmospheric models remain largely untested at 
non-solar metallicities, and there are no known benchmark brown dwarf companions 
to stars with significantly super- or sub-solar chemical compositions ([Fe/H] > -1-0.3 or 
[Fe/H] < -0.3). The latest-type ultracool subdwarf companion known is the d/sdM9 
benchmark HD 114762B ([Fe/H]= -0.7); atmospheric models do areasonably good job 
of reproducing the medium-resolution (A/AA ~ 3800) near-infrared spectral features 
of this object, but fits to the low resolution {A/AA ~ 150) near-infrared spectrum are 
unreliable (Bowler et al. 2009) . 

Although atmospheric models are not yet grounded by brown dwarfs with known 
metallicities, trends in the models have provided qualitative indications of deviations 
from solar metallicity for a growing number of L and T dwarfs with peculiar spectra. 
The optical spectra of peculiar L dwarfs are marked most notably by enhanced metal- 
hydride and metal-oxide bands compared to normal L dwarfs of the same spectral type 
(Figure[3]l. These variations are likely caused by subsolar metallicities and possibly sup- 
pressed condensate formation. A reduced metallicity also increases collision-induced 
absorption by H2 (CIA H2), resulting in bluer NIR colors for a given optical spectral 
type. Cloud properties also influence the NIR colors of L dwarfs and there is no clear 
way to distinguis h clouds from a niild d eviation from solar metallicity from NIR colors 



or spectra alone (IBurgasser et al. 2008b. Among the ~ 20 known objects that make up 



this class of "blue L dwarfs" ( Kirkpatrick et al. 201Clf) . which is distinct from L subd 
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Figure 3. Optical spectra of L subdwarfs (red). The most notable differences in 
metal-poor L dwarfs compared to ordinary field objects are a stronger CaH absorp- 
tion band at 6800 A and stronger TiO absorption bands at 7100 A and 8 400 A . 
From top to bottom the optical spectra originate from 'Burgasser et al.' ('20091), 
iBurgasser et al. (2007), Lodieuetal. (2010),"Cushing et al. (2009), Bowler etal] 
(120101). and [B urgasser et al.l (l2003l) . Comparison spectra (black) are from 
iKirkpatrick etal.. (1999); from top to bottom they are 2MASS 1146-H2230 (L3), 
2MASS 1155-1-2307 (L4), DENIS-P J1228.2-1547 (L5), DENIS-P J1228.2-1547 
(L5), 2MASS 0850+1057 (L6), and DENIS-P J0205.4-1159 (L7). The spectra are 
normalized between 7900 A and 8200 A and are offset by a constant. 
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Figure 4. The r - / vs. g-r color-color digram for stars in the SDSS spectroscopic 
catalog. The thick colored lines show the mean loci for the 4 metallicity classes of 
M dwarfs (dM:red; sdM:green; esdM:blue; and usdM :purple). 



warfs, two benchmark blue L dwarfs provide important clues about the nature of the 
spectral pecuharities. The blue L dwarf 2MASS J17 114559+4028578 (L4.5) orbits a 
solar-metallicity star (Radiaan et al. 2008) and the blue L dwarf SDSS J141624.08+ 134826.7 

(d/sdL6, iBowler et al., 2010. : Schmidt et al..,2010a) has a peculiar T7.5 compan i on with 

spect ral features indicative of being mildly metal-poor ( Bumingham et al.]|2010l : lBurgasser et al 
l2010h : this is the first evidence that blue L dwarfs may span a range of metallicities. 
For T d warfs, gravity an d metallicity both affect the /T-band flux by influencing CIA 
H2 (e.g.. iLiu et al.ll2007h . Metallicity (and to a lesser extent gravity) al so affects the 
F-ban d flux, offering a way to distinguish between these parameters (e.g.. Legget t et al.l 
l2007h . Ongoing sensitive all-sky surveys like WISE and Pan-STARRS are expected to 
greatly increase the census of non-solar metallicity isolated and benchmark L and T 
dwarfs, enabling rigorous testing of atmospheric models and an empirical calibration 
of spectral classification schemes. 



5. The Colors and sub-Classes of Subdwarfs in SDSS 



Low-mass stars with very low metallicities, typical of the Galactic thick disk and halo 
population, have a spectral energy distribution that is significantly different from the 
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Figure 5. Luminosity vs. temperature diagram for a series of isochrones with 
masses between 0.09 and 1 Mq and ages between 10^ and 10^ years with [Fe/H]=0. 
The red lines have Solar C and O abundances, while the blue and green lines are 
enhanced in O and C respectively. 



more metal-rich disk stars. The reason lies in the reduced absorption from metal oxide 
bands, in particular TiO. M dwarfs are classified in four so-called "metallicity classes" 
based on the relative strengths of their TiO bands: from the metal-rich dwarf M dwarfs 
(dM), to subdwarfs (sdM), extreme subdwarfs (esdM), and the ve r y me tal-poor ultra- 
subdwarfs (usdM). The classi fication follows the system of iGizisI ( 19971) recently up- 
graded by lLepine et al.l (l2007h . The seque nce usdM^esdM— »sdM— »dM is belie ved to 
form a sequence of increasing metallicity dGizis & Reidl[T997l : IWoolf et alJl200i) . with 
[Fe/H]w-0.5 for sdM, [Fe/H]w-1.0 for esdM, and [Fe/H]<-1.5 for usdM, although the 
metallicity calibration remain relatively uncertain to this date. 

A recent search of the Sloan Digital Sky Survey (SDSS) spectroscopic database 
(Lepine et al. in preparation) has produced over 7,600 M subdwarfs. Their color dis- 
tribution reveals significant differences with the metallicity class. Figure |4] shows the 
r - i color as a function of g - r. The dots show a typical distribution for nearby field 
stars, displaying the well-known "elbow" with a sti^ong inflection point at ^ - r = 1.4, 
r - i = 0.6. The thick colored lines show the mean loci for the dM (red), sdM (green), 
esdM (blue), and usdM (purple). There is a clear segregation as a function of the 
metallicity class, which reflects the strong effect that the TiO bands have on the spec- 
tral energy distribution. Ultrasubdwarfs simply extend the linear relationship between 
g - r and r - /, as one would expect from a blackbody. As the metallicity increases, 
the "elbow" becomes increasingly pronounced. This happens because the r-band gets 
increasingly depressed in the more metal rich stars, as the TiO opacity increases. This 
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strong dependence of color on metallicity opens the possibility of estimating metallici- 
ties in low-mass stars based on broadband photometry alone. As it turns out, even dM 
show a significant scatter in g - r which could be entirely explained by differences in 
metallicity. Should this be confirmed, this would provide a formidable tool for quick 
and easy metallicity estimates of low-mass stars. A proper calibration of the g - r and 
r - / color terms as a function of metallicity should be a priority. 



6. Metallicity and Stellar Evolution Models 

"Metallicity" loosely describes the heavy element content of a star or stellar population. 
Metallicity and [Fe/H] are often used interchangeably, with the implicit assumption that 
the other heavy elements scale with Fe as they do in the Sun. If they don't, then it's 
important to understand how changing a given element alters the spectrurt i, hence the 
opacity, hence the effective temperature scale of the star ( Potter et alj|2007 ). 



After H and He, the t wo most abundant e lements in the sun by mass or number 



fraction are C and O (e.g. jAsplund et al.ll2009l) . When C or O is enhanced relative to 



solar at fixed [Fe/H] the most dramatic effect appears in the molecular- opacities. Figure 
[5] shows a series of isochrones with masses between 0.09 and 1 Mq and ages between 
10^ and 10^ years with [Fe/H]=0. As Figure [5] indicates, on the one hand, enhancing 
C actually makes the lowest mass stars hotter while, on the other hand, enhancing O 
makes them cooler. This behavior can be understood in terms of the contribution of 
water molecules to the opacity. When analyzing the physical properties of low mass 
stars with effective temperatures below about 4,000K it is important to consider that 
non-solar abundance ratios can skew the results. 



7. Metallicity and Atmosphere Models 

One of the primary tools for measuring metallicities of low-mass objects is the compar- 
ison of data to synthetic spectra (c reated using atmo spheric models). Calculating such 
models is common practic e (e.g.. ISordo et al.l [20101) using modern atmosphere codes 
(iHauschildt & Baronl ll999h. For low-mass objects in general, it is crucial to account 



for dust formation in the atmosphere. This dust i: ormation needs to be treated as a 
microphysical growth and destruction process (e.g. iHelling et al.ll2008l) . Furthermore, 



molecules, both as an opacity source and as material that affects the equation of state, 
need to be accounted for with accurate input data such as formation constants or related 
quantities and line lists or equivalent opacity data. The situation is further complicated 
when c alculating models specifically for low metallicity objects, such as the models of 



Witte et al.i (,2009) . The decreasing metal content does not change or even simplify any 



of the main physical processes, but, in contrast, adds another dimension of parameter 
space. In particular , dust keeps formi ng in significant amounts down to metallicities of 
about [Fe/H] =-4.0 dWitte et al.ll2009l) . 

Recently, it has become possible to apply s ynthetic spectra to obs ervations of L 
subdwarfs and to attempt to measure metallicities ( !Burgassg r et al. 2009). However, the 
quality of the fits and the derived metallicities still vary (I witte et al.ii201Q) depending 
on the quality of the implemented physics (see also Fig. O. However, it is worth noting 
that the derived me tallicities of the sdL class do not need to be the same of the sdM class 
as derived by e.g. iGizisI d 19971) or lSchweitzerl (fT999i) . Measuring metallicities directly 
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Figure 6. The sdL4 dwarf 2MASS 1626+3925 (black, dotted: lBurgasseii l2004 
and a fit with the DRIFT 2009 models (blue, solid; Witte et al. 2009). The compaii- 
son with the DRIFT 2010 models (red, solid; Witte et al. 2010) is a comparison with 
the same model parameters showing the differences in model details. Both models 
have Teif = 2100K, log(g)=5.0 and a metalHcity [Fe/H]= -1.5. 



at a resolution of 0.1 dex or higher has not been attempted yet since the molecular 
background Unes add too much uncertainty. 



8. Conclusions 

During the first half of the last century, spectroscopic observations and radiative transfer 
theory began to unlock the composition of stars. Determining the metallicity of stars 
has been very important to a wide range of astronomical investigations, from planetary 
to cosmological scales. Yet, despite the progress made for most of the main sequence, 
measuring the metallicity of the Galaxy's most populous members, M dwarfs, remains 
a daunting task. 

At the start of this century, astronomers are beginning to unlock the metal content 
of these stars. However, there is much work to be done by both observers and theorists. 
Observationally, there are promising new results suggesting that IR observations may 
be important for estimating metallicities. This is strengthened by the relative agreement 
between models and spectra in this regime. However, these methods need further test- 
ing (with M dwarf binaries or clusters). In the optical bandpass, the relative metallicity 
classes described in Section 5 display a clear separation in photometric colors. This 
will be crucial for estimating the metal content of these stars in the next generation of 
surveys, which will be largely photometric. Yet, these classes have not been rigorously 
tied to an absolute metallicity scale. Once this occurs, the chemical composition of M 
dwarfs will be a powerful tool for studying the Galaxy and identifying the most likely 
exoplanet hosts. Identifying new benchmarks, for both M dwarfs and brown dwarfs, 
will be crucial in calibrating optical observations. 
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On the theoretical front, new line lists, opacity calculations and the inclusion of 
dust grains have resulted in better agreement with observations. The effects of carbon 
and oxygen abundance differences can be modeled and explain observations of the 
lower main sequence of globular clusters. As computational power and techniques 
advance, these models will grow in sophistication and should offer a more realistic 
picture of the important physics within these stars. 

Unlocking the metallicity of M dwarfs will profoundly benefit the astronomical 
community in a variety of ways, such as identifying exoplanet hosts and studying chem- 
ical evolution. The work presented here represents the first steps in solving this prob- 
lem. 
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